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A theory is presented for the binding of small molecules such as surfactants to semiflexible poly- 
r--s , mers. The persistence length is assumed to be large compared to the monomer size but much smaller 

f^ ■ than the total chain length. Such polymers {e.g., DNA) represent an intermediate case between flex- 

^SJ ' ible polymers and stiff, rod-like ones, whose association with small molecules was previously studied. 

The chains are not flexible enough to actively participate in the self-assembly, yet their fluctuations 
induce long-range attractive interactions between bound molecules. In cases where the binding 
significantly affects the local chain stiffness, those interactions lead to a very sharp, cooperative 
association. This scenario is of relevance to the association of DNA with surfactants and compact 
proteins such as RecA. External tension exerted on the chain is found to significantly modify the 
binding by suppressing the fluctuation-induced interaction. 
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'"' I. INTRODUCTION 

Ti ■ 

rj , Aqueous solutions containing polymers and small associating molecules such as folded proteins and amphiphiles 

Q ■ (surfactants) are commonly found in biological systems and industrial applications. As a result, extensive efforts have 

O , been devoted in the past few decades to the study of polymer-surfactant interactions |Q,||. In addition, there has 

been growing interest in the interactions between DNA macromolecules and surfactants, lipids and short polyamines 

r/-\ ' |B-10|. These interactions are relevant to various biochemical applications such as DNA extraction and purification 

J> 1^ 10| and genetic delivery systems ||ll|]. Association of folded proteins {e.g., RecA) with DNA plays a key role in 

C^ ' genetic regulatory mechanisms. Structural details of this association have been studied in recent experiments |12| , p^ . 

^«0 , Recently, we have presented a general theory for the self-assembly in aqueous solutions of polymers and smaller 

J^ ■ associating molecules |1J,0. Two different scenarios emerge, depending on the flexibility of the polymer. If the 

^•^ ' polymer is flexible enough, it actively participates in the self-assembly, resulting in mixed aggregates jointly formed 

^y. • by the two species. The polymer conformation changes considerably upon self-assembly but remains extended on a 

qn^ \ global scale, as the chain undergoes only partial collapse |l4|-[l6|. On the other hand, if the polymer is stiff, partial 

*>••; . collapse is inhibited. 

C^ ■ The criterion determining the 'flexible' vs. 'stiff' scenarios concerns the polymer statistics on a mesoscopic length 



o 



scale characterizing correlations in the solution (usually a few nanometers). It was found p 15| that the flexible (stiff) 
scenario holds if the exponent v, relating the number of monomers N to the spatial size R they occupy, R ~ N^ , is 

' ^ . smaller (larger) than 2/d on that length scale {d being the dimensionality). This distinction is analogous to the one 
made in the critical behavior of certain disordered systems ]1^7|l3l — if the critical exponent i/ of a system satisfies 
u < 2/d, the critical behavior would be smeared by impurities (in analogy to the partial collapse), whereas if i^ > 2/d, 
the critical point remains intact. Indeed, neutral flexible polymers in three dimensions, having v ~ 3/5 < 2/3, are 

.5^ found by scattering experiments to associate with surfactants in the form of a 'chain of wrapped aggregates' ||l9| , p0| . 

^^ , On the other hand, stiff DNA molecules, having i/ = 1 on the relevant length scale, are found to either remain 
\—{ ' unperturbed by surfactant binding jg,^] , or undergo a discontinuous coil-to-globule transition Q , provided the chain 

5t . is nruch longer than the persistence length. 

In previous publications |14| , |l5| we concentrated on the flexible case and the corresponding partial collapse, where 
the polymer degrees of freedom play an important role. In the opposite extreme limit of stiff, rod-like molecules, 
the conformational degrees of freedom of the polymer can be neglected and the chain may be regarded as a linear 
'binding substrate'. Models for stiff polymers, inspired by the Zimm-Bragg theory |M|, treat the bound molecules as 
a one-dimensional lattice-gas (or Ising) system with nearest-neighbor interactions p^]. They have been widely used 
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to fit experimental binding isotherms for polyelectrolytes and oppositely charged surfactants [^. Recently, more 
detailed electrostatic models have been proposed for the interaction between rod-like polyelectrolytes and oppositely 
charged surfactants p4|j25|l. In addition, a theoretical work focusing on the specific binding of proteins to DNA has 
been recently presented |[26[ , treating a pair of bound proteins as geometrically constraining inclusions on the DNA 
chain. 

In the current work we address the intermediate case of semiflexible polymers. The polymer we consider is stiff 
in the sense defined above, i.e., its persistence length, Zp, exceeds several nanometers and, hence, the polymer is 
characterized by z^ = 1 > 2/3 on that length scale. The total chain length, however, is considered to be much 
larger than /p, and the entire polymer cannot be regarded, therefore, as a single rigid rod. This case corresponds, in 
particular, to experiments on long DNA molecules p-[lO[, whose persistence length is typically very large (of order 50 
nm), but much smaller than the total chain length (which is usually larger than a micron) p^. We argue that such 
an intermediate system may, in certain cases, be governed by different physics. Although the polymer is too stiff to 
change conformation and actively participate in the self-assembly, its degrees of freedom induce attractive correlations 
between bound molecules. Those fluctuation-induced correlations are weak but have a long spatial range (of order Zp) 
and, hence, may strongly affect the binding thermodynamics. 

The model is presented in Sec. O. Bound molecules are assumed to modify the local features of polymer confor- 
mation, e.g., change its local stiffness. In the limit of weak coupling, our model reduces to the Kac-Baker model 



8|-|30[, which is solvable exactly. This limit is discussed in Sec. HI. Although turning out to be of limited interest 



in practice, the weak-coupling limit provides insight into the mechanism of association, and helps us justify further 



approximations. Section IV presents a mean-field calculation for an arbitrary strength of coupling. This analysis leads 
to our main conclusions, and in Sec. M it is extended to polymers under external tension. The results are summarized 



in Sec. VI, where we also discuss several relevant experiments involving DNA and point at future directions. 



II. THE MODEL 

Small molecules bound to stiff polymers are commonly modeled as a one-dimensional lattice gas (or Ising system) 
^2| . Each monomer serves as a binding site, which can either accommodate a small molecule or be empty, and 
the surrounding dilute solution is considered merely as a bulk reservoir of small molecules. In the current work we 
stay at the level of a one-dimensional model, assuming that the polymer is still quite (yet not infinitely) stiff, i.e., 
the persistence length is much larger than the monomer size. In addition, a dilute polymer limit is assumed, where 
inter-chain effects can be neglected. We focus on the effect of introducing the polymer degrees of freedom and, hence, 
seek a simple meaningful coupling between the polymer and the bound 'lattice gas'. 

A polymer configuration is defined by a set of vectors, {un}n=i...N, specifying the lengths and orientations of the 
A'^ monomers. In addition, each monomer serves as a binding site which can be either empty ((^„ = 0) or occupied by 
a small molecule ((^„ = 1). A configuration of the entire system is defined, therefore, by specifying {u„, (Pn}n=i...N- 

Since the polymer is assumed to be locally stiff, a natural choice would be to couple (/?„ with the square of the local 
chain curvature, (p„(u„^i — u„)^, thus modifying the local chain stiffness. However, in the usual Kratky-Porod worm- 
like-chain model of semiflexible polymers [Q, chain segments are taken as rigid rods of fixed length (|u„| = const), 
and each squared-curvature term contains only one degree of freedom {e.g., the angle 9n between u„ and u„+i). 
Consequently, this coupling, ip„ cos0„, would leave {(fn} uncorrelated, leading merely to a trivial shift in the chemical 
potential of bound molecules [p3] . One option to proceed is to consider higher-order extensions of the worm- like-chain 
Hamiltonian, involving three consecutive monomers. This will introduce correlations between bound molecules at 
different sites. 

We take, however, a simpler route and modify the worm-like-chain model by allowing the monomer length to 



fluctuate. This modification was originally presented by Harris and Hearst |33|, using a single global constraint for 
the average chain length. The modified model was shown to successfully reproduce the results of the Kratky-Porod 
model as far as thermodynamic averages {e.g., correlation functions, radius of gyration) were concerned. It was less 
successful, however, in recovering more detailed statistics of the worm- like chain {e.g., distribution function, form 
factor), particularly in the limit of large stiffness. The Harris-Hearst model was later refined by Lagowski et al. |3J] 
and Ha and Thirumalai [p5lB6|, replacing the single global constraint by a set of local constraints for the average 
segment lengths. This further modification was shown to be equivalent to a stationary-phase approximation for the 
chain partition function, yielding reliable results for average quantities, as well as more detailed statistics [B5|. We 
note that a similar approach was used in a recent model of semiflexible polymer collapse [|37| . It should be borne in 
mind that, despite its success in the past, the constraint relaxation remains essentially an uncontrolled approximation. 
In the current work we restrict ourselves to thermodynamic averages, such as monomer-monomer correlations and 
free energies, for which the modifled model with a single global contraint can be trusted. 



Thus, the rigid constraints of the original Kratky-Porod model, u^ = 1, are relaxed into thermodynamic-average 
ones, (u^) = 1, where the mean-square monomer size is taken hereafter as the unit length. Using the modified model 
for the chain, each ((5„(u„_|_i — u„)^ term involves two consecutive monomers (and not merely the angle between them), 
leading to a meaningful coupling between binding and polymer conformation. 

The partition function of the combined system of polymer and bound molecules is written, therefore, as 



Z = 



N 

Tr / Y\ du„ exp(-H) 



JV-l Af N 

W = T^p ^ (1 + e(/J„)(u„+i - U„)^ + ^ Kul - fi^ifn- (1) 

n—1 n—1 n—1 

In Eq. (|l|) Zp is the persistence length of the bare chain, characterizing its intrinsic stiffness. It is assumed to be much 
larger than the monomer size, Zp 3> 1. The coupling is introduced through the stiffness term, assuming that a bound 
molecule modifies the local stiffness by a fraction e > — 1, which may be either negative or positive but cannot change 
the positive sign of the overall stiffness term p8[. The second term contains a set of multipliers, A„, to be chosen so 
that the constraints (w^) = 1 are satisfied. However, replacement of the entire set {A„} by a single multiplier A can be 
shown to yield a non-extensive correction [pq], which becomes negligible in the limit N ^ oo. Hence, we use hereafter 
a single multiplier, A. Finally, the system is assumed to be in contact with a reservoir of solute molecules. The last 
term in Eq. p) accounts for this contact along with any other factors which couple linearly to the degree of binding. 
Typically, /i contains the chemical potential of the solute reservoir and the direct energy of solute molecule-monomer 
binding. (All energies in this work are expressed in units of the thermal energy, k-QT.) Note that we have not included 
in Eq. (|]) any direct short-range {e.g., nearest-neighbor) interactions between bound molecules. Thus, all interactions 
in the model arise from the coupling to the polymer degrees of freedom. Short-range interactions between bound 
molecules do exist in physical systems. Yet, in the limit of Zp ^ 1 and |e| ^ 1, which is of interest to the current work, 
such direct interactions have a minor effect on binding, as is demonstrated in the following sections. Hence, we omit 
them for the sake of brevity. 

As a reference, let us start with the previously studied partition function of the bare polymer [pq] , 

Zp = / J|du„exp[--/p^(u„+i-u„)2- A^w^]. (2) 

n n n 

It is a Gaussian integral which can be calculated either by transforming it to Fourier space and integrating, or by 
analogy to the path integral of a three-dimensional quantum oscillator p9|] . The result in the limit N ^ oo and for 

(3) 



Zp > 1 is 




/ , \ 3/2 / 


\ 


^r^CJ -(- 


- V3^/'pj ■ 


The multiplier A can now be determined according to 




1 dlogZp 2\ 1 





(4) 

where (• • •)p denotes a thermal average over the bare chain statistics {i.e., using Zp). The corresponding free energy 
per monomer (in the ensemble of constrained u„) is 

1 3 3 

/p = -— logZp- A = -log/p + — +const. (5) 

Various correlations in the bare chain can be calculated. The pair correlation between segment vectors along the 
chain sequence is 

(u„ • u„)p = e-l™-"l/'^ (6) 

which explains why the parameter Ip has been defined as the persistence length. Two higher-order pair correlations 
are calculated as well: 

gi = ((u„+i - u„)2)p ^ + o{lp^) 
tp 

52(m, n) EE ((u™+i - u„0'(u„+i - u.^)\ ~ gj = J^e-2|'"-"l/'p + 0{l;'), (7) 

and will be of use in the next section, where we re-examine the coupled system. 



III. WEAK COUPLING 

Let us return to the full partition function m), which can be equivalently written as 

Z = Zp Tr exp(^^(^„) /exp[ 1- ^ (/?„(u„+i - u„)^] ) • (8) 

n \ n 



First we consider the weak-coupling limit, |e| <C 1, where the partition function (g|) can be treated by a cumulant 
expansion. In this limit the model becomes analogous to the exactly solvable Kac-Baker model ES-BG] , and we show 
that identical results are derived from a simple mean-field calculation. We then use this observation to justify a 
mean- field calculation for an arbitrary value of e. 

A cumulant expansion of Eq. (g) to 2nd order in e leads to 



Z ~ Zd Tr exp 



S^nE \ x-^ 1 /3Le^ 



M I~-9i ) X!'^"+ 9 ( ^ ) ^92{'m,n)Lp„,Lp. 



(9) 



where the correlations gi and (72 were defined in Eq. (|7|). Substituting expressions (|7|), the partition function is 
decoupled into a polymer contribution and an effective contribution from the bound solute molecules, 

Z ~ ZpZs = Zp Tr exp(-7Ys) 

Ws = ^ X! ^"iriVm'Pn - fi^V,i, (10) 



where 



2 , 



V = _^p-2|m-n|/ip 

3e 3e2 

The introduction of the polymer degrees of freedom and their coupling to the binding ones have led to two effects, as 
compared to previous lattice-gas theories. First, there is a shift in the chemical potential, /i — > /i. This is equivalent to 
an effective change in the affinity between the small molecules and the chain. As expected, if binding strengthens the 
local stiffness of the chain (e > 0), the affinity is reduced {i.e., the isotherm is shifted to higher chemical potentials), 
whereas if it weakens the stiffness (e < 0), the shift is to lower /i. The second, more interesting effect is that bound 
molecules experience an attractive potential, Vmm along the chain. The amplitude of this effective interaction is small 
(~ e^/lp), but its range is large — of order Ip. When Ip is increased there are two opposing consequences — the 
interaction amplitude diminishes, while the interaction range is extended. The overall effect on the thermodynamics 
of binding, therefore, has to be checked in detail. 

A. Analogy with the Kac-Baker Model 

The effective Hamiltonian of the bound solute. Tig, is a lattice-gas version of the Kac-Baker model p8f-p0[, which 
is exactly solvable. Moreover, the procedure relevant to our semiflexible polymer, i.e., increasing Ip while keeping 
1 -^ Ip -^ N, is precisely the one studied in detail by Kac and Baker. Their results, as applied to our binding 
problem, can be summarized as follows. For any finite Ip, the bound molecules are always in a disordered state 
along the polymer chain, as in any one-dimensional system with finite-range interactions. Consequently, the binding 
isotherm, i.e., the binding degree (p = (fn) as function of /i (see, e.g.. Fig. 0a), is a continuous curve. However, in the 
limit Zp — > cxD, taken after the infinite-chain limit N —^ oo, there is a critical value of coupling above which the binding 
exhibits a discontinuous (Ist-order) transition. According to Baker's rigorous calculation [Q, the critical value of the 
potential amplitude multiplied by Ip (equal, in our case, to 3e^/2) is 4, i.e.. 



±1.63. (12) 



Note that the symmetry with respect to the sign of e is merely an artificial consequence of our 2nd-order expansion, 
Eq. (^ . In general, the results should not be the same if the stiffness is weakened (e < 0) or strengthened (e > 0) , as 
is demonstrated in Sec. pv| . 

The negative critical value in Eq. (|lj), e~ ~ —1.63, lies outside the range of validity of the original polymer binding 
model, e > — 1 [cf. Eq. (fi^)]. The positive value, e+ ~ 1.63, does not satisfy the assumption of weak coupling, 
|e| <C 1, which have led to the analogy with the Kac-Baker model in the first place. Thus, the sharp binding isotherms 
obtained from the Kac-Baker model for |e| > ec do not apply, strictly speaking, for our polymer binding problem. The 
weak-coupling calculation does demonstrate, however, how fluctuations in polymer conformation induce long-range 
attraction between bound molecules. This basic feature is expected to remain when one considers stronger coupling, 
|e| > 1, and the resulting many-body terms omitted in Eq. (^). This is further discussed in the following sections. 

Finally, the polymers we consider have a large but finite Ip. For example, the persistence length of a DNA macro- 
molecule is typically of order 50-100 nm, whereas the length of a single base pair is 0.34 nm. Hence, Ip is of order 10^ 
(in units of monomer length) . It is worth checking to what extent the sharpness of binding in the Kac-Baker model 
for |e| > Ec is affected by finite Ip. For this purpose, let us define a cooperativity parameter ior the binding, measuring 
the maximum slope of the binding isotherm. 






(13) 



This parameter is equivalent to the zero magnetic field susceptibility in the analogous spin system, and is commonly 
measured from the slope of binding isotherms obtained in potentiometric experiments [0,0. It has been defined in 
Eq. (g_3|) so as to yield zero for vanishing interaction (e = 0) and diverge at a critical point. (In the current weak- 
coupling limit, the maximum slope is obtained for {(p) — 1/2.) Given Ip and e, the cooperativity is numerically 
calculated using Kac's exact solution |22j2£|, as is explained in the Appendix. Figure || presents the results for Ip — 10 
and 50. For Ip = 50 the binding becomes highly cooperative for |e| > ec- For even larger values of Ip ~ 10^ (relevant, 
e.g., to DNA) the binding will be hardly distinguishable from that of an infinite Ip. 

B. Mean-Field Calculation 

In fact, the results of the Kac-Baker model in the limit N -^ oo,lp -^ cxd, while keeping Ip < N, can be also 
obtained from a simple mean-field calculation |2q,Bn| . The heuristic argument for this agreement is the following: as 
Ip is increased, the range of interaction is extended and each bound molecule interacts with an increasing number of 
neighbors. As a result, the averaging assumption underlying the mean-field approximation is justified, and becomes 
exact when the range of interaction is taken to infinity. The correspondence between infinite-range models and mean 
field was rigorously proved by Lebowitz and Penrose for a more general class of potentials ma . 

Indeed, employing a mean-field approximation for the potential ( |ll|) in the limit of very large Ip, 

mn ^ \ mn / 

where ip is an average, uniform binding degree, we are led to the following mean-field free energy per monomer: 

f ^ fp + fs^ fp + 'f iog If + {l-ip)\og{l- ip)- —ip'^ - flip, for^p^oo. (14) 

It is easily verified that the critical point of this free energy is e^ = 8/3, in agreement with the rigorous result, Eq. (p2). 
The cooperativity parameter can be calculated as well from Eq. (14), yielding 

^ = 4(g2 _ g2) ' fo"- Ip^oo. (15) 

This expression shows the usual critical behavior obtained from mean-field theories, C ^ \e ~ ed"''' with 7 = 1. The 
dependence of C on e according to Eq. (|lj) is shown by the solid line in Fig. ||. The curves obtained from Kac's 
solution approach it, as expected, when ^p is increased. Recall that expressions ( p^ ) and ( p^ correspond to the 
original problem of bound molecules only in the limit of small e. 



IV. STRONG COUPLING 

The interesting part of our theory requires |e| J^ 1 and thus hmits the interest in the analogy to the Kac-Baker 
modeL Nevertheless, based on the heuristic argument given above, it is reasonable to assume that, in the limit 
Ip^ 1, the mean- field approximation is good for larger values of |e| as well [Q. The preceding section, discussing the 
Kac-Baker model in the weak-coupling limit, may be regarded, therefore, as a justification for using the mean-field 
approximation for one-dimensional models with large Ip and |e| ^ 1. Applying a mean-field approximation to the 
binding degrees of freedom ipn in our starting point, Eq. dif), the tracing over u„ can be done exactly. The resulting 
free energy is composed of the polymer free energy, /p, evaluated with an effective persistence length, Ip — > /p(H- eip), 
and the entropy of mixing for ip, 

f = fp\i^^i^{i+ev) +V^ogip+{l^ip) log(l -if)- flip. (16) 

Using Eq. (^, we obtain 

/ = <^ log (^ + (1 - ^) log(l -'/')+ 2 log[^p(l + (-ip)] + ^^ .^_^ ^ - ^^V■ (17) 

For e <C 1 and /p ^ 1 this expression reduces, as expected, to our previous result for the weak-coupling limit, Eq. (M). 
In the limit Ip ^ 1 the critical points of the free energy ( p7| ) are 

e- - ^ (2 - x/To) ^ -0.775, e+ = ^ (2 + Vio) ~ 3.44, (18) 

both of which lie within our general range of validity, e > — 1 . (Note the loss of symmetry with respect to the sign 



of e, which was a consequence of the weak-coupling approximation in Sec. HI.) The corresponding critical chemical 
potentials are 

f^t - ^H^ - l°g(^c^ + 1) - ±1-67. (19) 

The binding isotherm, ip — ipi^^, as derived from Eq. (O), satisfies 



ip 3e 

1~ ^ 2(1 + £</?) ' 



y. = log ^^— + :^7^^3, Ip » 1. (20) 



Figure Ba shows three binding isotherms for three different values of e below and above the critical point. The 
corresponding binding cooperativity is 

C ^(1 + ^)' 1 Z»l (21) 

^'3(2 + .)2(.-.c-)(e+-.) 4' 'P»^- ('^^ 

As in Eq. (iq), this expression exhibits the usual mean-field critical behavior, C ^ |e — ed"''' with 7=1. The 
dependence of C on e is plotted in Fig. ||b. 

Finally, the binding phase diagram arising from Eq. ( p7| ) in the limit /p ::^ 1 is depicted in Fig. 0. At the lower 
limit of model validity, e -^ —1, the spinodal approaches a finite value, /isp = log(2/3) — 5/2 ~ —2.91, whereas the 
binodal diverges. Indeed, for e ^ — 1 the free energy ( p7| ) tends to —00 for (/? = 1, regardless of the value of /i, and the 
binodal is thus obtained at /i ^ —00. In this respect, the limit e = — 1 for the bound molecules is similar to the limit 
of zero temperature — the induced interaction is so strong that the molecules condense for any value of the chemical 
potential. Note that in this special limit, e -^ — 1, </3 ^ 1, the effective stiffness, lp(\ -I- eip), becomes vanishingl y sm all. 
This limit cannot be accurately treated within the continuum form of the semifiexible polymer Hamiltonian |p8| . 

Equations (18)^(|2l|) and the phase diagrams in Fig. summarize the results obtained so far. They indicate that in 



cases of semifiexible polymers, where binding of small molecules significantly affects local chain features, the binding 
should be a very sharp process. For finite Ip the slope of the binding isotherm is finite, i.e., the binding is always 
continuous, yet for Ip ^ 10^ like in DNA, the behavior will be practically indistinguishable from a discontinuous phase 
transition. 

It should be borne in mind that the sharp binding, obtained despite the one-dimensionality of the model, relies on 
the long range of the induced interaction. A direct short-range interaction between bound molecules could not produce 
a similar effect. Hence, such a short-range interaction (e.g., a nearest-neighbor interaction), which was omitted in 
Eq. (nl) for the sake of brevity, does not have an important effect on the binding in the domain of interest, i.e., Zp S> 1 
and lei > 1. 



V. CHAINS UNDER TENSION 

In addition, we consider binding to semiflexible chains which are subject to external tension. This scenario is 
relevant to recent single-molecule manipulation experiments |T2,n3|. Since the tension suppresses chain fluctuations, 
it is expected to have a significant effect on the fluctuation-induced mechanism presented in the preceding sections. 

In order to incorporate the external tension into our model, a term is to be added to the chain Hamiltonian [cf. 

Eq. (I)] ii, 

N 



Z^ Tr / TT du„ exp(-H - Ht 

N 
Ht = -t • ^ U„, 



(22) 



ra=l 



where Ti. has been defined in Eq. (|l|), and t is the exerted tension (in units of ksT divided by monomer length). 

As in Se c. O , we begin with the previously studied problem of a bare semiflexible chain, yet it is now a chain under 
tension |36| , |43|] . The additional tension term has not changed the Gaussian form of the polymer part of Z. It can be 
calculated, therefore, in a similar way to that of Sec. ||, yielding 

Z^/"' = Zy^exp(iV4A), (23) 

where Zp is the tensionless polymer partition function given in Eq. (j^). The equation for the multiplier A is, in this 



case, 



^pA 



1/2 



4A 



which reduces to Eq. m) for t = 0. The resulting polymer free energy is 



/pt = ^^oglp 



1/2 



4A-^' 



(24) 



(25) 



where A — X{lp,t) is the solution to Eq. 
For Ipt -^ 1, the solution for A is 



and the free energy becomes 



4^ 



4?, 



-(iptf + oiiptY 



fp,c,fp~jf + 0{iy), t^l/lp 



(26) 



where /p is the tensionless free energy given in Eq. (H). This is the elastic regime, where the energy is quadratic 
{i.e., the relative chain extension is linear) in tension |3q|44[| . Since we assume a large persistence length, this regime 
corresponds to very weak tension, t <C l//p <C 1. In the opposite limit, Ipt '3> 1, the solution to Eq. (04) becomes 



A 



and the corresponding free energy is 



/. 



pt 



2^oglp 



3t 



2lr 



2lpt 



1/2 



1/2 



oiipt)-' 



*+Ur +^(^pn, t->i/ip 



(27) 



In this regime the chain extension changes like the inverse square root of tension pq,Ea . 

Let us turn now to the effect of tension on the system of polymer and bound molecules, Eq. (p^). As in Sec. IV, we 
employ the mean-field approximation, valid for /p — > cxd. The resulting free energy is the same as Eq. (p7|), but with 
/pt instead of /p. 



/= fpt\ i^^i^^i+,^) + V^og if + {1 - if) log{l - ip) - ^iip. (28) 

Due to the additional degree of freedom, namely tension, the binding phase diagrams of Fig. y become three- 
dimensional. In particular, the critical points ej become critical lines, ej(i). (Note that t is an external field 
coupled to {u„} rather than {ifin}, and, hence, it does not destroy the critical behavior.) The 'condensation' of 
bound molecules in our model results from attraction induced by polymer fluctuations. By suppressing fluctuations, 
the tension should weaken the attraction and shift the critical coupling to higher values, i.e., increase the positive 
critical point, e+, and decrease the negative one, e~. Using Eqs. (|24|), ( P5| ) and (P8|), the critical lines, e^(i), can be 
calculated. The results are shown in Fig. ^. 

Before getting into the detailed effect of tension, we address the question whether the critical behavior can survive 
any strength of tension. In this respect there is an essential difference between stiffness-strengthening binding (e > 0) 
and stiffness- weakening one (e < 0). In the former case, since the value of e is unbound, there exists ejt(i) for any 
value of t, such that the binding is a sharp transition for e > e^(i). In other words, the critical line e^(i) exists for 



any < t < cx3. Indeed, substituting e ^ oo in Eq. (|2S| ) while using Eq. (27), we find that the free energy always 
describes a sharp transition, regardless of the value of t. On the other hand, in the latter case of stiffness-weakening 
binding, there is a lower bound for e, e > — 1, where the validity of the entire approach breaks (see previous section). 



Substituting e = — I in Eqs. (Eq) and (27), we find that a critical point exists only for t < t* , where 

^ (S3- 7721) ~ 0.410. (29) 

Thus, the critical line e~{t) terminates at the point (i*,e* — ~1), beyond which a sharp binding transition cannot 
be attained. This situation is similar to a case where the critical temperature Tc coincides with T = (e.g., in a 
one-dimensional Ising model), and the system is disordered at all temperatures T > 0. 

Several regimes are found as function of tension. For very weak tension, t < 1/lp, the leading-order term which 
couples binding and tension is found from Eqs. (Eq) and (Eq) to scale like Ipt^etp, i.e., it is only linear in 93. Hence, 
to leading order in Ipt there is no effect on the critical point. Although the tension influences chain fluctuations {e.g., 
causing the chain to extend linearly with i), it is too weak to affect the fluctuation- induced interactions between 
bound molecules. The next-order term scales like l^f^i^ + ^'P)^j leading to a very small shift of ^ l^f^ in the critical 
point (see also Fig. |^). 



For t > 1/lp, the leading-order term in the free energy, according to Eqs. (27) and (p8|), is ~ (i//p)"'^/^(l -I- eip) ^1"^ . 
Here two regimes should be distinguished. For intermediate tension, 1/Zp < t < Ip, the critical line scales like (i/Zp)^/^, 
reflecting a more significant, yet still weak effect of tension. Although the chain conformation is significantly stretched 
by tension in this regime, the induced interaction between bound molecules is not strongly affected. However, for 
t > Ip, the tension term in the free energy [~ (i//p)^/^(l + eip)^^^^] becomes dominant, leading to a linear dependence 
of the critical point on tension, e^ ~ t/lp. 

The above analysis for the dependence of the critical coupling on tension is summarized in the following expression: 

l^f^ t < l/lp 

let{t)-etm ^ { (/p/i)i/2 l//p<t<;p (30) 

Ip/t t > Ip, relevant only to e+. 

The various regimes are also clearly seen in Fig. ^. Note that for the large values of Ip considered in this theory the 
intermediate tension region, 1/Zp < t < Ip, is very wide. 

VI. DISCUSSION AND CONCLUSIONS 

We have considered binding of small molecules to isolated semiflexible polymer chains, where the persistence length 
Ip is much larger than the monomer size but still smaller than the total chain length N. We have demonstrated that 
in such systems polymer fluctuations induce attraction between bound molecules. The long range of this interaction 
(of the same order as the persistence length) can lead to strong effects on the binding process. In particular, if 
bound molecules signiflcantly affect local features of the chain, e.g., weaken or strengthen the stiffness by a factor 
of about 5 (e < e~ or e > e+), then the binding is predicted to be extremely cooperative, occurring as a transition 
for a sharply deflned solute concentration. This is an unusual, yet practical example for a one-dimensional system 
exhibiting a sharp transition due to long-range interactions. The results of the model should apply, in particular, to 
the association of DNA with smaller molecules such as surfactants and compact proteins. 



Subjecting the polymer to external tension has been studied as well. By suppressing the fluctuation-induced 
interaction, the applied tension may strongly affect the binding. The effect is significant for sufficiently strong tension 
of order t ~ Ip. [For DNA this implies t ^ 10^fcBT'/(10A) ~ 10^ pN.] In cases where binding weakens the chain 
stiffness, such a high tension should make the sharp binding transition disappear altogether (i.e., regardless of the 
strength of coupling or temperature). In cases where binding strengthens the chain stiffness, a tension of i ^ /p 
significantly shifts e^ to higher values. It is worth mentioning that tension-induced pairwise interaction between 
specifically bound proteins on a DNA chain was studied in a previous work pq | . 

The interaction of DNA with oppositely charged cationic surfactants has been thoroughly studied by potentiometric 
techniques |^,y and fluorescence microscopy |^,||. Isotherms measured by potentiometry reveal a very cooperative, 
albeit continuous binding. Fluorescence microscopy convincingly demonstrated, however, that the binding to a sin- 
gle DNA molecule has a discrete nature resembling a Ist-order phase transition. It is usually accompanied by a 
coil-to-globule collapse of the DNA chain (which lies outside the scope of the current theory). The smoothness of 
potentiometric isotherms was shown to arise from averaging over an ensemble of DNA molecules, coexisting in bound 
and unbound states g . Similar results were obtained for the association of DNA with spermidine |0] . The microscopic 
origin of the observed cooperativity (or even discontinuous transition) has not been clarified. It is usually fitted to a 
phenomenological parameter describing strong interaction between nearest- neighboring bound molecules 122] . On the 
other hand, it is reasonable to expect that oppositely charged surfactants bound to DNA chains significantly modify 
the chain stiffness (probably weakening it). Thus, our model demonstrates that the strong cooperativity observed in 
experiments can be well accounted for by weak, yet long-range interactions induced by polymer fluctuations. 

Recently, the kinetics of non-specific binding of RecA proteins to DNA has been studied by single-molecule manip- 
ulation ||l2|,|l3[. RecA is a bacterial protein involved in DNA recombination and known to cause significant changes 
in the local structure of the double strand upon binding El] . It was found to increase the DNA stiffness by a large 
factor, estimated around 10 in one study ||l^ and above 4 in another |]1^. This corresponds to a large, positive e 
in our model. A very cooperative nucleation-and-growth kinetics was observed, as expected from the current model. 
Moreover, in certain situations it was possible to achieve a smaller increase of stiffness by binding of RecA. This led, 
correspondingly, to a less cooperative process |13| ]. Yet probably the most compelling evidence is that the binding 
cooperativity was shown to be sensitive to external tension of order 10-100 pN. It was consequently deduced that 
DNA conformational fluctuations play a key role in RecA binding ||l^ , in accord with the model. 

The current work is restricted to one-dimensional interactions along the chain sequence, assuming that the polymer 
is locally stiff and obeys the worm-like-chain description. Apart from changing local properties of the polymer, an 
important feature not treated by the model is that bound molecules may also modify volume interactions between 
the monomers, thus affecting the three-dimensional conformation of the polymer. For example, binding of oppositely 
charged surfactants to a DNA molecule locally neutralizes the DNA charge. This should lead, indeed, to a modified 
stiffness, but also to a reduced 2nd virial coefficient, which may drive a coil-to-globule collapse g. The collapse can 
be also driven by fluctuations in the concentration of ions adjacent to the chain, as has been demonstrated by recent 
theoretical studies ]37|j47|] . 

In order to check the theory presented in this work more experiments are required, focusing, in particular, on 
the effect of persistence length and tension on binding. The fluorescence microscopy techniques, which have been 
successfully used for DNA-surfactant association, may be applied to chains under tension or flow, thus examining 
the role of fluctuations. It would be interesting to study a system consisting of a semiflexible polymer and bound 
molecules in computer simulations, and thereby check the applicability of our mean-field approximation. An important 
extension of the model, as mentioned above, would be to introduce volume interactions and obtain binding-induced 
collapse as observed in experiments. 
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APPENDIX: NUMERICAL DETAILS 

The aim of the numerical scheme is to calculate the results of the Kac-Baker model for finite l^, which are presented 
in Fig. |l|. Using Kac's solution ^8 29|, the partition function of bound solute molecules, Eqs. ([lO|)-([ll|), is expressed 
in the limit A'^ ^ oo as 



Zs = const X e^ 



J 



where cq is the largest eigenvalue of the following 'transfer kernel': 



K{x, y) = [(1 + e^-3e/2+VJ:.^(^ ^ gM-3e/2+Vjy^ji/2 



exp 



y'-^' 



(y 



— o ^/'pr: 



2(i_e-4/;p) 



(Al) 



(A2) 



where J = 3e^/2Zp, and x,y E (—00,00) are real variables. 

We define a vector, {xi] = {{2i — M)d}i=o,,,M, where M is an even integer and d a real number, and use it to 
discretize the kernel K{x,y) into a transfer matrix. 



^j — ^ y-^i^ -^j ) • 



In addition, we define the diagonal matrix 



XiO'ij . 



(A3) 



(A4) 



Given Ip, e and /i, the transfer matrix Kij is diagonalized and its largest eigenvalue, eo, is found. 

The binding degree, 1^9, can be calculated in two ways. The first is by calculating the variation of log eo with respect 
to /x, 



The second way is by using the equation 



Lp^dloge^/dn. 



if - ioo/(5\/j), 



(A5) 



(A6) 



where B = coth(l/^p), and A is the matrix A transformed to the basis where K is diagonal [J48| . 

The cooperativity parameter, C, as defined in Eq. (|l3h, is found by calculating the variation of <~p with respect to 
/x around the point ip = 1/2. The value /i = M1/2 which gives Lp — 1/2 is analytically found by transforming the 
lattice-gas partition function, Eqs. (pX|)-(pl|), into an Ising one {ipn -^ s„ = 2(^„ — 1), and requiring that the 'magnetic 
field' coefficient should vanish. The result is 



Ail/2 = 3e/2- JB/2. 



(A7) 



For each calculation (i.e., for each set of Zp, e and /i) the discretization parameters, M and d, were tuned until the 
result became insensitive to further refinement to six significant figures. In addition, the two methods for calculating 
ip were used and verified to yield identical results to six figures. All algebraic manipulations were performed using 
Mathematica. 
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FIG. 1. Binding cooperativity as function of e according to the Kac-Baker model, plotted on a semi-logarithmic scale. The 
dashed and dash-dotted curves are results of numerical calculations for /p = 10 and 50, respectively. The solid curves show 
analytic results for Ip 
lines) . 



cxD as obtained by a mean-field calculation [Eq. (Ilq)]. The critical points are at e^ = ±y/8/3 (dotted 
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FIG. 2. (a) Binding isotherms obtained from the mean-field theory [Eq. (12(1)] for three different values of e: e — 1 (solid 
line), £ = 3 (dashed) and e — 4 (dotted), the latter being beyond the critical point ej ~ 3.44. The chemical potential /i 
is given in units of ksT. (b) Binding cooperativity as function of e according to the mean-field calculation [Eq. (Ell)]. The 
cooperativity diverges at the two critical points e^ ~ 2(2 4z \/T0)/3 (dotted lines), beyond which binding isotherms exhibit a 
Ist-order transition [see dotted curve in (a)]. 
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given by Eqs. (hq) and (191). The chemical potential /i is 
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FIG. 4. Effect of tension on the critical behavior of binding, (a) Critical coupling, etr(i) < 0, as function of tension. Two 
regimes are found: for t < l/'p, |er(i)| increases like t^; for t > l/'p, it increases like i}'"^ . The critical line terminates at the 
point (t* ~ 0.410/p,ec ~ ^1), beyond which a sharp binding transition becomes unattainable, (b) Critical coupling, e^(t) > 0, 
as function of tension. Apart from the two regimes of (a) there is a third one for t ^ Zp, where ej|'(i) increases linearly with t. 
The value taken for /p in the numerical calculation is 100. The tension t is given in units of fceT divided by monomer length. 
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